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Abstract The Swift/XRT light curve of the ultraluminous X-ray (ULX) source NGC 5408 
X-l was re-analyzed with two new numerical approaches, Weighted Wavelet Z-transform 
(WWZ) and CLEANest, that are different from previous studies. Both techniques detected 
a prominent periodicity with a time scale ofll5.5±1.5 days, in excellent agreement with 
the detection of the same periodicity first reported by Strohmayer (2009). Monte Carlo 
simulation was employed to test the statisiticak confidence of the 115.5-day periodicity, 
yielding a statistical significance of > 99.98% (or > 3.8er). The robust detection of the 
115.5-day quasi-periodic oscillations (QPOs), if it is due to the orbital motion of the 
binary, would infer a mass of a few thousand M Q for the central black hole, implying an 
intermediate-mass black hole in NGC 5408 X-l. 
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1 INTRODUCTION 

The presence of supermassive black holes (SMBHs, 10 6 — 10 9 Mq) in the hearts of active galactic nuclei 
(AGNs), and the stellar-mass black holes (BH) (~40 M Q ) as the evolutionary remnants of massive 
stars has been widely accepted and also confirmed by extensive observations with various techniques. 
However the existence of intermediate-mass black holes (IMBHs, 10 2 — 10 4 M Q ) remains a puzzle. 
IMBHs are thought as the seeds of the SMBHs, the masses of which were subsequently built up either 
through the accretion of gas, or through hierarchical black hole mergers (reviewed bv lVolonteri 20101 ). 
Therefore discovery of IMBHs is of vital importance for understanding the dynamics of stellar clusters, 
formation of SMBHs and production of gravitational waves dMiller & Colbert 20041 ). 

The detection of a population of ultraluminous X-ray sources (ULXs, X-ray luminosity exceed- 
ing 10 39 erg s _1 ) invokes an interpretation that they might be accretion-powered objects which har- 



bour IMBHs (e.g. Colbert & Mushotzky 1999 1. However, the high X-ray luminosity may also result 



from relativistically beamed emission. In this regard, the central object is not necessary an IMBH 



(King et al. 2001 1 , or the high X-ray luminosity may represent a super-Eddington accretion of stellar- 



mass black holes (Begelman 2002). Alternatively some regular ULXs might represent an extension of 
low-mass X-ray binaries (LMXBs; ILiu et al. 20061 ). Determination of the mass of ULXs is crucial to 
distinguish between these rival interpretations of the physical nature of ULXs. The canonical methods 
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of measuring dynamical mass that rely on direct measurement of the binary motion of the companion 
are often frustrated because the companion is usually too dim to be detected. On the other hand, indirect 
methods that are based on the fit of X-ray spectral energy distribution (e.g., |Makishima et al. 2000t and 
X-ray timing observations (e.g. JKaaret et al. 20061) also place constraints on the black hole mass . 

NGC 5408 X-l is one of the few ULXs showing variable X-ray flux (ISoria et al. 20041). The 
source shows QPOs around 10-40 mHz observed with the XMM-Newton ( jStrohmayer et al. 2007| 



Dheeraj & Strohmayer2012i. Recently Strohmayer (2009) claimed to discover QPOs in NGC 5408 



X-l with P = 115. 5 ± 4 days from the Swift/XRT monitoring data. However Kaaret & Feng (2009) 
analyzed the same data set (the only difference is that Strohmayer did exposure map corrections in the 
light curve) but did not detect significant periodicity, although they indeed found a clue of a 115-day 
periodicity. In contrast with Galactic X-ray binaries whose characteristic variability period is of the or- 
der of days, it requires a rather long time span to search for periodicity of ~ 1 00 days in NGC 5408 X- 1 . 
Except for M82 X-l ( |Kaaret & Feng 2007] ), there are few other ULXs that have been monitored for such 
a long period. The QPO of NGC 5408 X-l is crucial for better constraining the physical nature of this 
ULX: if the 1 15.5-day periodicity is due to the orbital motion of the binary, that would provide an evi- 
dence for an IMBH in NGC 5408 X-l ( jStrohmayer & Mushotzky 2009|>. Both groups mainly relied on 
the periodicity analysis on the Lomb-Scargle (LS) periodogram dLomb 1976l|Scargle 1981) . However 
the statisitical significance of the LS-technique-detected periodicities is still a controversial issue. The 
LS periodogram introduces a traditional confidence limit based on an exponential distribution e~ z (z 
is the highest peak), that is actually not accurate enough for random signals. Therefore additional time 
series analysis techniques are necessary to verify the detected periodicity. 

In this Letter, we shall re-analyze the light curve of NGC 5408 X-l to check the 115.5-day peri- 
odicity by using two other numerical techniques, the Weighted Wavelet Z-transform (WWZ) and the 
CLEANest methods which are more robust than the LS periodogram. The statistical significance of 
detected periodicities will also be evaluated with the Monte Carlo simulations. 



2 THE ALGORITHM 



Traditional methods which are used for searching for and identifying periodic fluctuations in time se- 
ries are based on Fourier transform. A modified Fourier-based periodogram method is called Date- 
Compensated Discrete Fourier Transform (DCDFT: Fe rraz-Mello 198 H i which models the observed 
data as a linear combination of three trial functions of constant l(t), cos (tut), sin(ujt). Astronomical 
data are usually not evenly sampled. However Fourier transform of an unevenly spaced time series may 
often introduce a myriad of complications, resulting in altering the peak frequency slightly and chang- 
ing the amplitude greatly, and even generating strong fake peaks. In order to identifying true periodic 
signals out of a time series containing multiple periodicities, an updated algorithm named CLEANest 
JFoster 1995b was developed. The CLEANest technique successively subtracts the strongest peaks in the 
original power spectrum detected by the DCDFT, until the remaining strongest peak in the residual spec- 
trum is statistically not significant. These processes stop when the residual spectrum is considerably low 
and flat. The advantage of the CLEANest technique is that the significant frequencies involved would be 
refined at each iteration of the CLEANest procedure, thus the fake peaks induced by the sidelobe effect 
of the primary peaks are greatly reduced. 

Compared to Fourier analysis method which represents only frequency analysis, the wavelet trans- 
form is remarkable for its localization properperty in both time and freqeuncy domains. The wavelet 
transform decomposes the signals into a combination of wavelet functions which can be scaled and 
shifted in time and frequency. The (quasi-)periodic behaviors of the signals can be identified through 
visual observations of the wavelet power as a function of translation (time) and dilation (frequency), 
and the characteristic frequency (or time scale) is determined by looking for peaks in the time-averaged 
wavelet power (Grossm ann et al. 19891 ). For unevenly-sampled data, the fluctuations of the local num- 
ber density may induce spurious high-frequency spikes and also result in fake time evolution behav- 
ior of the characteristic periods. In order to solve these problems, Foster (1996) proposed a rescaled 
wavelet method called Weighted Wavelet Z-transform (WWZ). Same with the canonical wavelet trans- 
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Fig. 1 WWZ periodogram of NGC 5408 X-l (a) and the statistical histograms of Monte- 
Carlo simulations of the 210-day periodicity (b) and of the 1 15.5-day periodicity (c). The red 
dashed lines in (b) and (c) denote the peak powers of the periodic components in Fig. 1-a. To 
the right of the boundary line, the WWZ power of the simulated light curves is higher than 
the observed value; the ratio of the accumulated count of these light curves to the total test 
number (i.e., 5000) corresponds to a probability of the periodicity induced by the pure noise. 



forms, the WWZ is also a kind of function projection which projects the signals onto three trial func- 
tions: a constant function (fi(t) = l(t), <f2(t) = cos(ui(t — t)), 933(4) = sin(u(t — r)), where u 
is the dilation parameter, also called scale factor, and r is the time shift.. But the difference is that 
the WWZ adopts statistical weights u) a = exp[—cuj 2 (t a — t) 2 ] (o=1, 2, 3) to the projection on the 
trial functions, where the constant c determines how rapidly the wavelet decays. The WWZ power is 
defined as ^-statistics, WWZ = (N eff - 3)/2(V x - V y ) (Foster 1996), where N eff is the effec- 
tive number of the data which represents the statistical data density, and V x and V y are the weighted 
variation of the data and the weighted variation of the model function, respectively. In doing so, the 
WWZ follows the F-distribution with N e f / — 3 and 2 degrees of freedom. The CLEANest and WWZ 
technique have been widely used in studies of variable stars (e.g., |Kiss etal. 19991 Templeton 2004 1. 
Recently they are also successfully applied in X-ray timing analysis of AGNs and X-ray binaries (e.g., 
|Espaillat et al. 2008|ILachowicz et al. 20061 1. 



3 RESULTS 

The periodicity analysis results derived from the WWZ are presented in Figure Q] The WWZ power is 
shown as a function of both the observing time (x-axis) and the testing period (y-axis) (Fig.QJa). The 
peaks in the WWZ power indicate the strength and duration of a characteristic periodicity. The majority 
of the WWZ power is dominated by the most prominent component whose peak marks a characteristic 
period of 1 15.5 days, and this periodic component spreads across the whole time span with a scattering 
of ±1.5 days. The amplitude of the 1 15.5-day component is 14.2, significantly higher than the secondary 
peak of 6.0, which corresponds to a charasteristic period of ^210 days. The ridge line of the contours 
of the 210-day periodic component shows a large variation from 184 to 226 days. 

Figure 2 shows the power spectra made from the Fourier transform analysis: the DCDFT (Fig. 2-a) 
and the CLEANest (Fig. 2-b). In the DCDFT spectrum (Fig. 2-a), there are two peaks higher than the 
half of the maximal amplitude. The primary peak is located at a characteristic frequency of 0.00866 
day -1 , corresponding to a period of 115.5 days. The secondary peak falls in a relatively broad range 
from 194 to 236 days; the highest peak denotes a period of ~210 days. 
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The periodograms derived from the WWZ and DCDFT show excellent consistency, and they 
are in agreement with the previous results based on the Lomb-Scargle periodogram analysis 



( Strohmayer 2009 1 as well: a distinct primary peak at a period of 1 15.5 days, and a scattered secondary 
at a period of ~210 days. Note Strohmayer (2009) also detected the secondary peak at P2~210 days in 
his Figure 2, although the author did not remark on the secondary. The harmonic relation between two 
characteristic periods (frequency ratio ^2:1) and the scattering appearance of the secondary are remi- 
niscent that the 210-day periodicity might be an artifact due to the sidelobe of the primary. Therefore we 
further run the CLEANest algorithm to check the reliability of the secondary periodic component. After 
subtracting the strongest 1 15.5-day peak, the residual CLEANest spectrum shown in Fig. 2-b becomes 
rather flat, and does not show any significant peak any more. That gives evidence that the 210-day peak 
detected in the WWZ and DCDFT spectrum is probably a fake signal which is contaminated by the 
nearby primary 1 15.5-day periodic component. 

In order to further constrain the nature of the periodicities detected by the WWZ, we performed 
statisitical confidence tests with the Monte Carlo simulations dLinnell Nem ec & Nemec 1985). First we 
created a series of 5000 light curves following Possion noise distribution to represent the random ob- 
servational errors. The simulated light curves had the same number of the sampling (113), temporal 
baseline (485 days), mean sampling rate (~14 days) and variance as the observed light curve. The ob- 
served data were randomized in this way. We then calculated the WWZ periodograms for each of the 
simulated light curves using the same parameters as we did for the observed light curve. Next in the fre- 
quency channel correponding to the 1 15.5-day periodicity, we searched for the periodograms in which 
the peak power exceeds 14.2, i.e., the peak power of the 1 15.5-day periodicity, and counted the number 
of the light curves. This accumulated number represents the occurrence that the 115.5-day periodicity 
can be generated by pure noise. The ratio of this occurrence to the total testing number denotes a prob- 
ability of fake detection of the periodicity due to the noise in the observed light curve. The 210-day 
periodicity was tested following the same procedure. The histogram of the Monte Carlo simulations are 
displayed in FigureQJb (210-day) andQJc (1 15.5-day). In Fig. 1-c, most of the power of the simulated 
light curves is lower than 8.0, there are only two exceptions higher than 10.0 but lower than 14.2 (the 
peak power derived from the observed data), placing a lower limit of the confidence level of 99.98% 
(or > 3.8cr). Compared with the observed WWZ periodogram in Fig. 1-a, the simulations provide a 
unambiguous evidence for the detection of the 1 15.5-day periodicity across the whole observing period. 
In Fig. 1-b, there are 13 simulated light curves that may generate peak powers of the 210-day periodicity 
higher than that from the observed data, setting a upper limit of the statisitical confidence of 99.74% 
(< 3a). 
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Fig. 2 The DCDFT (Top) 
and CLEANest (Bottom) pe- 
riodograms of NGC 5408 X- 
1. The red solid line indicates 
the 1 15. 5 -day periodic compo- 
nent. The red dashed line indi- 
cates the fake 210-day period. 
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4 DISCUSSION 

In the present work, two numerical techniques, the WWZ and the CLEANest, were used to search for 
the periodicities in the Swift/XRT light curve of NGC 5408 X- 1 . Both WWZ and DCDFT periodograms 
detected a distinct component across the whole time span correponding to a 1 15.5-day period, in good 
agreement with previous works ( |Strohmayer 2009[ ). Monte Carlo simulation was employed to test the 
statisitical significance of the 115.5-day periodicity. The simulations indicate a high confidence level 
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of >99.98%, excluding the possibility that the periodicity is generated by the random noise. All the 
evidence points to a conclusion that the X-ray light curve of NGC 5248 X-l intrinsically displays QPOs 
with a time scale of 1 15.5 days. 

A straightforward interpretation of the 115.5-day periodicity is the orbital period of the binary 
system, resembling with the QPOs detected in Galactic BH binaries such as XTE J 1550-564 and GRS 
1915+105 (|Strohmayer 2009| and references therein). The difference between NGC 5408 X-l and the 



Galactic BH binaries is that the time scale in NGC 5408 X-l is longer by a factor ^100, inferring a larger 
dynamic mass of the central BH. If the low-frequency QPO in NGC 5408 X-l is analogous to Galactic 
black hole binaries in physical nature, the mass of the NGC 5408 X-l black hole might be estimated 
from the Galactic BH mass and the scaling factor of the QPO frequencies. That gives an estimated mass 
of NGC 5408 X-l BH ranging from 1000 M Q to 9000 M Q ( |Strohmayer & Mushotzky 2009} . 

In addition to the orbital motion scenario, there are other arguments for this QPO in NGC 5408 
X-l. Foster et al. (2010) claimed that the 115.5-day QPO is probably a manifestion of jet precession 
resembling SS 433, and is thus superorbital in nature. Therefore NGC 5408 X-l might in fact be a 
large stellar-mass BH accreting with a super-Eddington rate to account for the high X-ray luminosity 
(IMidd leton et al. 201 11 1. Similarly, Soria (2007) proposed that the accretion properties of ULXs are dif- 
ferent from stellar-mass BHs, and low-frequency milli-Hz QPOs cannot be used as an indicator of the 
high black hole mass. In their model, the ULXs can be 50-100 M Q -mass BHs with high accretion rate. 

Continuing X-ray timing observations and energy spectra studies are necessary to discriminate be- 
tween the scenarios for explaining the 1 15.5-day periodicity. If further observations can prove that the 
1 15.5-day periodicity is indeed related to the orbital motion of the binary system, NGC 5408 X-l would 
be one of the few IMBH candidates in extragalaxies. 
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